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ABSTRACT 

Dust-obscured galaxies (DOGs) are a UV-faint, infrared-bright galaxy population that reside at 
z ~ 2 and are believed to be in a phase of dusty star-forming and AGN activity. We present far- 
infrared (far-IR) observations of a complete sample of DOGs in the 2 deg^ of the Cosmic Evolution 
Survey (COSMOS). The 3077 DOGs have (z) = 1.9 ± 0.3 and are selected from 24 ^m and r+ 
observations using a color cut of r+ — [24] > 7.5 (AB mag) and ^24 > 100 /^Jy. Based on the mid-IR 
spectral energy distributions, 47% are bump DOGs (star-formation dominated) and 10% arc power- 
law DOGs (AGN-dominatcd). We use SPIRE far-IR photometry from the Herschel Extragalactic 
Multi-tiered Survey (HerMES) to calculate the IR luminosity and characteristic dust temperature 
for the 1572 (51%) DOGs that are detected at 250 urn (> 3a). For the remaining 1505 (49%) that 
are undetected, we perform a median stacking analysis to probe fainter luminosities. Detected and 
undetected DOGs have average luminosities of (2.8 ± 0.4) x IO^^Lq and (0.77 ± 0.08) x lO^^Lg, and 
dust temperatures of (34 ± 7) K and (31 ± 3) K, respectively. DOGs contribute 30% to the 24 /im- 
selected galaxies' infrared luminosity function, calculated using far-infrared observations. If the IR 
luminosity is extrapolated using the 24 /im flux density alone, it is overestimated by a factor of 2 
on average. DOGs contribute 10 — 30% to the total star formation rate density of the Universe at 
z = 1.5 — 2.5, dominated by 250 /im detected and bump DOGs, compared to around 40% for all 24 jim 
galaxies above our flux limit. DOGs have a large scatter about the star-formation main sequence and 
their specific star-formation rates show that the observed phase of star-formation could be responsible 
for their observed stellar mass at z ~ 2. 
Subject headings: Galaxies: luminosity function — Galaxies: star formation — Infrared: galaxies 
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1. INTRODUCTION 

The bolomctric luminosities of Luminous Infrared 
Galaxies (LIRGs; Lm > lO^L©) and Uhra-LIRGs 
(ULIRGs; LiR > IO^^Lq) arc dominated by re- 
processed thermal dust emission, due to prodigious 
rates of star-formation activity, black hole accretion, 
or both. Locally, these sources are rare, although 
out to z ~ 1 they become more numerous and in- 
creasingly dominate the infrared lumi nosity function 
of galaxie s with i ncreasing reds h ift (iLe Floc'h et aD 



2005; Pe rez-Gonzale z ct al .. 120051: iCaouti et all 120071: 
Magnelh et al.ll2009t iRodighiero et al.ll2010l : lEales et al.l 
2010f ). (U)LIRGs are thought to trace a phase of in- 
tense star-formation activity, which is likely followed 
by, or partially concurrent with, an episode of vigorous 
black hole accretion. Upon the cessa tion of these phases 
each produce s an early- type galaxy (Gcnzcl ct al."2001t 
Farrah et all 120031 [Lonsdale et al. 2006; Veilleux et t^ 
20091) . 

Studies using the Multiband Im aging Photometer for 
Spitzer (MIPS; iRieke et all I2004D instrument on board 
the Spitzer Space Telescope (jWerner et al.l l2004| ) have 
identified high-redshift ULIRGs from their 24 /xm emis- 
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sion (lYan et al.ll2004l: iHouck et al.ll200 5HWee dnian et al 



2006; 



20081 



Fiore et al.l 



Fiore et alj 



200a iDev et alj 1 2008: F arrah et al 



20091 ) . iDevetaLi (|2008i ) exploited 



this technique in the Bootes field of the NOAO Deep 
Wide Field Survey (NDWFS) and presented a sample 
of dusty ULIRGs selected by the color cut R - [24] > 14 
(Vega magnitudes; 82^/8 r > 1000). This color selection 
sainples the rest-frame 7.7 fiui polycyclic aromatic hy- 
drocarbon (PAH) feature found in star-forming galaxies 
and causes the redshift distribution to have a biased 
average at z '-^ 2. Also, at 2 ^ 2, the DOG selection falls 
within range of the power-law component of AGN emis- 
sion in the mid-IR, which also identifies a population of 
active galactic nuclei (AGN). It is proposed that these 
dust-obscured galaxies (DOGs) are the latter stage of 



the sub-millimeter galaxy (SM Gs: iHuehes et al 



iSmailet al.lll997t iBarger et all 119981 iBlain et al 



1998; 



19991 



among others) phase where an AGN is triggered while 
star formation is still occurring, causing some du st to be 
heated to higher temperatures (IDev et al.l 20081) than in 
classic 850 fiui selected SMGs. iPope et alJ (|2008[) found 
that 30% of the SMGs are also DOGs, and of those 
SMG-DOGs, 30% are AGN-dominatcd (> 50% AGN 
contribution in mid-IR), consistent with this scenario. 
Using high resolution optical and near-IR (NIR) imaging 
from the Hubble Space Te lescope to invest i gate DOG 
morpholo gy, the s tudie s of iBussmann et al] (|2009f ) and 
JBussman n et al.l (j2011| ) found that the morphologies 
of bump (star-forming) DOGs, power-law (AGN domi- 
nated) DOGs, SMGs, and high redshift quiescent distant 
red galaxies (DRGs) are consistent with the picture that 
major merger-driven systems event ually all evolve into 
compact relaxed passive galaxies (iSpringel et al.l 12005 
and re ferences therein). Furthermore, iNaravanan et ah 
(|2010( ) used N-body and hydrodynamic simulations to 
model the temporal evolution of high redshift galaxies 
and found that at the peak of the merger-driven galaxies' 
star formation rate, a galaxy can both be identified as 
an SMG and a DOG. The same study also found that 
during the stages after final coalescence, merger-driven 
DOGs transition from being star-formation dominated 
to being AGN dominated. 



The launch of t he Herschel Space Observatory^^ 
(jPilbratt et al.l l20To( ) enables the direct observation of 
DOGs in the far-IR regime, instead of extrapolating from 
spec tral energy distr i bution (SEP) teinplates or stacking 
(e.g.lDev et ai.ll200l iPope et al.ll2008l ). lMelbourne et al.l 
()2012[) studied _ffersc/ieZ-detected DOGs with spectro- 
scopic redshifts and showed that DOGs classified by 
their mid-IR SEDs as cither bump (star-forming) or 
power-law (AGN-dominatcd) have 250 /^m/24 /im flux- 
density ratios that are con sistent with l o cal U LIRGs 
of the respective classes. iPenner et al.l (j2012| ) used 
Herschel data to show that DOGs' high rest-frame 
MIR/UV flux density ratios are due to varying amounts 
of UV dust obscuration, and speculated that it is caused 
by differing degrees of alignment between dust and stars, 
or simply by the differences in total dust content. 



^^ Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consortia 
and with important participation from NASA. 



The focus of this paper is to extend the study of DOGs 
in the far-IR to a complete sample in order to extract 
statistical conclusions about this galaxy population. 
We generate our DOG catalog using Subaru r+ band 
and MIPS 24 /im da ta from the Cosmolo gical Evolution 
Survey (COSMOS; IScoviUe et ail \2(M ) and combine 
it with multi-wavelength data in the far-IR from the 
Herschel Multi-tie red Extragalactic Survey^^ (HcrMES; 
lOliver et al.ll201^ . We calculate IR luminosities, star 
formation rate (SFR) and dust temperatures for all 
DOGs detected at 250 fiui and employ a stacking analy- 
sis to calculate the average properties of the undetected 
population and thus to probe fainter luminosities. For 
DOGs at z = 1.5 — 2.5, we generate a luminosity function 
and calculate the star-formation rate density at z ^ 2. 

This paper is organized as follows. Section |2| de- 
scribes the dataset and sample selection. The results 
and our analysis are presented in Section [31 We sum- 
marize our conclusions in Section |4l Unless specifi- 
cally stated, all magnitudes are reported in the AB 
system, where — 2.51ogjQS'y(/iJy) -I- 23.9 = AB mag, 
and assume a standard ACDM cosmology with Hq — 
70 km s"^Mpc"\ f^M = 0.3, and ^a = 0.7. 

2. DATA AND SAMPLE SELECTION 

2.1. Far- Infrared Data 

The 250, 350, and 500 fiui far-IR data were obtained 
using the Herschel -Spectral and Phot o metric Imaging 
Receiver (SPIRE; Gr iffin et all IMot iSwinvard etal] 
|2010[ ) as part of HerMES, with an area coverage that 
completely overlaps with the MIPS observations of the 
2deg2 COSMOS field. We use the first data release 
(DRl) of Her MES maps that were processed using the 
smap pipeline (jLevenson et al.ir2010[ ). The reduced maps 
reach 3cr point source depths of 8, 10, and 14 mJy, in 
the 250, 350; and 500 /im channels respectively, where a 
is the combined instrumental and confusion noise. For 
sources with 6*250 ^ '^'^i wc use the listed photome- 
try from the HerMES cross- identification catalog (XID). 
This catalog uses known positions of 24 /im sources as 
a prior, and estimates SPIRE fiuxes via linear inver- 
sion methods. Model selections are used to account for, 
and prevent overfitting, and to optimize the 24 //m in- 
put. The fitting meth od is outlined in more detail in 
iRoseboomet all ()2010f ). 

2.2. Optical and Mid-Infrared Data 

We use deep Subaru Suprime-Cam (|Komivama et aLl 
[2003!) aperture-corr ected r"*" photomet ry supplied by the 
COSMOS catalog (jCapak et al. l l2007l ). The 5a point- 
source depth for a 3" aperture is 26.8 mag. 

The mid-IR data are from Spitzer observations car- 
ried out by the CO SMOS Spitzer Survey (S-COSMOS; 
iSanders et al.l 120071) using the Infrared Array Camera 
(IRAC: iFazio et all I200I and MIPS. The IRAC 5cr 
depths at 3.6, 4.5, 5.6, and 8.0 /im for an aperture- 
corrected 1.9" aperture, are 0.50, 0.6, 3, and 5 /iJy, re- 
spectively. The MIPS 24.0 fiva ha point source depth is 
55 /iJy. 

We next generate a MIPS24/im-selected catalog that 
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combines the Subaru and Spitzer datasets, using a two- 
step cross-matching process within the 2 deg^ of the Sub- 
aru deep are a in order to f ind optical counterparts for 
each source (jFu et al.ll2010( ). Firstly, the 24 /im coordi- 
nates are matched to the closest IRAC detection within 
a 2" search radius, then the nearest optical counter- 
part is identified within 1" of the IRAC position. Fi- 
nally, sources near bright stars that were within the Sub- 
aru/optical and Spitzer/IRAC 3.6 fim coverage were re- 
moved from the catalog to avoid contamination. The 
final catalog has 48,474 sources and is > 90% complete 
above S24 = 80 //Jy. 

2.3. Sample Selection 

DOGs are selected in the standard manner, by identi- 
fying sources with r+ - [24] > 7.5 (AB mag; 824/ Sr+ > 
1000) and we require [24] < 18.90 mag (524 > 100/.tJy) 
due to the depth of the 24 /im data. Using these crite- 
ria 3077 of the 48474 (6%) COSMOS 24 fim sources are 
identified as DOGs (Fig.[l]). The mid-IR SED of each 
DOG is examined using IRAC photometry (> 5a) to 
classify whether a DOG contains a bump-like feature or 
resembles a power-law. For this study, a "bump" DOG is 
defined if it satisfies one of the following: 6*3.6 £ 5*4.5 > 
5*8.0; 5*4.5 < 55.8 > 5*8.0; or 63.6 < 54.5 > 55.8. Here 
5'[3.6,4.5,5.8,8.o] represent the flux densities in the 4 IRAC 
channels. Conversely, we label a DOG as "power-law" 
if it satisfies 63.6 < ^4.5 < 65.8 < 58.o- Previous stud- 
ies have interpreted sources that feature a bump in the 
mid-IR SED to be the stellar continuum peak at rest- 
frame 1.6 /im, tracing ste llar emission and likely star- 
forma tion dominated fe.g. lYan et al.|[2005t iSaiina et al.l 
|2007() . while a pow er-law is dominated by AGN con- 
tinuu m emission (e.g. lWeedman et 311120061 : iDonlev et al.l 
|2007| ). Bump DOGs compose 47% of our sample, while 
power-law DOGs are rarer, totaling 10%. The remain- 
ing 43% are not classified due to one of two possibilities: 
insufficient or low signal to noise IRAC data; or an SED 
shape that does not satisfy the above criteria. For the 
latter case, most of the sources are at z < 2 (median of 
z = 1.1), such that the rest-frame 1.6 /xm stellar contin- 
uum peak lies outside the wavelength range of the IRAC 
channels. 

2.4. Redshifts 

All redshifts used in this paper are from COSMOS. 
Spectroscopic redshifts are us ed when available (35 
sources, 1%; iLillv et all 120071 Kartaltepe et al., in 
prep), although virtually all of our DOG sample (2979 
sources, 97%) use photometric redshifts. The photo- 
metric redsh i fts ar e derived from 30 photometric bands 
(|Ilbert et al.l l2009f ). providing crAz/(i+z) = 0.02, for 
24 /im sources that lie at z = 1.5 — 3 and have the same 
r+ mag range as DOGs. The 61 DOGs that are X- 
ray detected use photom etric redshifts that a lso account 
for AGN flux variability (jSalvato et al.ll2009[) . Also, two 
sources do not have a redshift estimate and are excluded 
from our sample. We note that the sharp peak in the 
redshift distribution at z = 1.95 is due to rounding from 
the redshift values associated with the bin size used and 
no spatial correlation is observed. 

The redshift distribution of the final sample of 3075 
DOGs is shown in Fig. [2l with a mean of {z) = 1.9 ± 0.3. 
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Figure 1. r"''-[24] as a function of 24 /im flux. DOGs arc selected 
to have r+ - [24] > 7.5 AB mag and ^24 > 100 fiiy. DOGs with 
z = 1.5 — 2.5 are highlighted in red, while green arrows are lower 
limits for sources that were undetected in the r^-band. The 1901 
DOGs at 2 = 1.5 — 2.5 are not biased in r+-[24] when compared 
to other DOGs. 
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Figure 2. Photometric redshift distribution of DOGs in t he COS- 
MOS flcld. We show the DOG distribution from Bussrna nn et al] 
II2OI2I ) normalized to have equal peaks for comparison. The filled 
region highlights the range z = 1.5 — 2.5, considered for our analy- 
sis. We find (z) = 1.9 ± 0.3, assuming a Gaussian distribution, as 
shown in red. 

The sample of 90 D OGs in the Bootes field with spectro- 
scopic redshifts from lBussmann et al.l ()2012j ). normalized 
to have an equal peak with our sample, is also shown. 
The two samples have a consistent mea n z of 1.9 ± 0.02 
and 2 .1 ±0.5, for our sample and the iBussmann et al.l 
(|2012| ) sample, respectively. 

3. ANALYSES AND RESULTS 

3.1. Far-Infrared Spectral Energy Distributions 

Using the COSMOS redshifts and Herschel 250, 350 
and 500 /im photometry, we fit the far-IR SED and cal- 
culate the rest-frame IR luminosity (8 — 1000 /im) and 
characteristic dust temperature. We divide the DOGs 
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Figure 3. Example SED fitting for a detected DOG {left panel) 
and an undetected DOG {right panel). The black curve shows the 
best fitting template to the SPIRE data points (black circles) and 
the gray curves show CEOl templates that provide acceptable fits 
consistent with the error bars. The red curve shows the best fit 
modified blackbody, which we use to calculate the dust tempera- 
ture. 



into two subsamples based on 250 ^m detections: a DOG 
is considered -ffersc/ieZ-detected if it satisfies S'250 > 3cr250 
(where (T250 is the total uncertainty due to the instrumen- 
tal and confusion noise), and undetected otherwise. Of 
our DOG sample, 51% are thus 77ersc/ie/-detected. To 
calculate the characteristic dust temperature, for each of 
these we use the available SPIRE flux densities to fit a 
modified blackbody of the form 



S^ ex B„{Tdust)i^' 



fi 



(1) 



where u is frequency, /3 is the dust e missivity, fixed to 
the typical value of 1.5 (jDrainel[2003D . Tdust is the dust 
temperature and Bi, is the Planck function, defined as 



B,= 



2hi 



1 



ahf/k-BT^u 



1 



(2) 



Here h is Planck's constant, c is speed of light, and ^b is 
Boltzniann's constant. The temperature we calculate is 
insensitive to and consistent with the reported error bars 
from varying /3 slightly. All 250 /xm-detected DOGs are 
also detected in one other band can be fitted with this 
function. 

We derive estimates of the IR luminosity by fitting the 
available SPIRE data to the SED template library of 
ICharv fc Elba3 (|2001f ) (hereafter CEOl). The template 
with the minimum x^ is chosen for the best fit. The un- 
certainty in IR luminosity is derived by first producing 
1000 mock catalogs for each source that assume a Gaus- 
sian distribution centered around the listed SPIRE flux 
density, with a dispersion equal to the average SPIRE 
flux density error. The IR luminosity per source is re- 
calculated 1000 times and the standard deviation of the 
IR luminosity distribution is the error in our calculation. 
Examples of the SED template and modiflcd blackbody 
fitting are shown in Fig. [31 

The IR luminosity (8 — 1000 um) is converted to star 
formation rate using (jKennicuttlll998[ ) 



SFR(M0yr- 



1.72 X 10-^"LiR(io) 



(3) 



which assumes a Salpeter initial mass function (IMF). 
We note that in our study we assume that UV emis- 
sion will provide negligible cont ribution to the total star- 
formation rate, as validated bv lPenner et al.l ()2012l) . 
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Figure 4. Example median stacking results of Herschel- 
undetected DOGs at 250, 350, and 500 /xm for undetected DOGs 
at 2 = 1.75 to 2.00 



To measure the average flux density of the undetected 
DOGs, we bin the sources in redshift and for each bin 
stack on the SPIRE residual maps. These maps are 
generated by performing a blind extraction and PSF- 
subtraction to prevent contamination of individually de- 
tected sources. We use the public l y ava ilable idl stack- 
ing library from iBethermin et al.l ( 20101 ) to perform the 
stacking^'^. Each stacked image was converted from the 
native Jy beam" to Jy pixeP and aperture photome- 
try with an aperture size equal to 22", 30", and 42" for 



250, 350 and 500 ^m respectively, is performed to calcu- 
late the flux of the stacked images. These aperture flux 
densities are consistent with those measured in the cen- 
tral pixel when the stacked map is in units of Jy beam~^ . 

The observed stacked flux densities are corrected for 
the boosting from clustering bias by dividing by fac- 
tors of 1.07, 1.10 and 1.20 at 250, 350 and 500 /im, re- 
spectively. The appropriate correction factors vary with 
clustering strength and are thu s population dependent . 
These values were calculated bv IBethermin et aLl (|2012[ ) 
for 24 ^m sources and are valid for DOGs because the ob- 
served correlation le ngths, rp (a proxy fo r clustering am- 
phtude), for DOGs (jBrodwin ey il."2008') and the parent 
population of 24 ^tm sources (jMa glio cch etti et al.l 120081 : 
iStarikova et al.ir2012f ) are consistent. Errors in the pho- 
tometry are calculated from bootstrapping the sources 
to be stacked. For each redshift bin, the clustering- 
corrected SPIRE flux densities of undetected DOGs are 
set to equal the median stacked flux densities and the IR 
luminosity and dust temperature are calculated using the 
same method as for the Herschel-detected DOGs. The 
(clustering-corrected) stacked fluxes and errors, and the 
resulting average infrared luminosities and dust temper- 
atures are shown in Table [T] 

We note that the stacks of 250 fim images are on av- 
erage a factor of 2 deeper than our 3cr detection limit. 
In Fig. in we show an example of the median stacked 
images for 250, 350 and 500 /im from left to right at 
z = 1.75—2.00 and an example SED using stacked SPIRE 
flux densities for an undetected DOG at z = 1.88 is 
shown in the right panel of Fig. |3l Each image stack 
is large enough to provide a good estimate for the back- 
ground noise. 

Fig. shows IR luminosities of the /fersc/ie/-detcctcd 
DOGs as a function of redshift. The average IR lumi- 
nosity for detected and undetected DOGs is (2.8 ±0.3) x 
10^^ Lq and (6.0 ± 1.0) x 10" Lq, respectively. LIRGs 



23 -jj^p YBL librarv from iBethermin et al. (2010?) is available at 
|http://www.ias.u-psud.fr/irgalaxics/downloads.php| 
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qa 
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qa 


N'' 




rpc 

dust 




(mjy) 


(mJy) 


(mJy) 




(XlOl^Lg) 


(K) 


< 1.50 


4.2 ±0.5 


3.0 ±0.4 


2.1 ±0.4 


354 


0.16 ±0.13 


25.1 ±5.5 


1.50-1.75 


3.5 ±0.5 


2.6 ±0.5 


1.7 ±0.4 


218 


0.37 ±0.02 


34.6 ±0.9 


1.75-2.00 


3.4 ±0.5 


2.9 ±0.4 


1.6 ±0.4 


406 


0.52 ± 0.07 


37.2 ±0.9 


2.00-2.25 


4.3 ±0.6 


3.9 ±0.6 


1.7 ±0.5 


237 


0.68 ±0.06 


40.0 ±1.0 


2.25-2.50 


4.3 ±1.1 


3.7 ±0.9 


2.7 ±0.7 


82 


1.00 ±0.21 


37.8 ±0.8 


>2.5 


4.9 ±0.7 


4.0 ±0.6 


2.9 ±0.6 


185 


1.71 ±2.8 


44.3 ±6.3 



Notes 

° Listed flux densities arc from median stacking 

*" Number of sources per redshift bin. 

'^ Average and standard deviation per bin. 

(10" < Lir(Lo) < 10^2) comprise 15% of Herschel- 
detected DOGs and 75% for Herschel-undetected DOGs. 
ULIRGs (10^2 < Lir(Lo) < 1013) make up 78% of 
the /fersc/ie/-detected and 15% for 7/ersc/ie^-undetected 
DOGs. Hyper-luminous infrared galaxies (HLIRGs (> 
10^3 ^0) are the rarest, totaling 2% for ffersc/ie/-detected 
DOGs and none for ifersc/ieZ-undetected DOGs. Al- 
though we note that there is additional uncertainty in 
the fractional contributions of the undetected sources, 
due to the use of stacked average fluxes, which mini- 
mizes the contribution from extreme sources. Herschel- 
detected power-law, or AGN-dominated DOGs, have on 
average Lir = (4.5 ± 0.5) x lO^^^ Lq, making them more 
IR-luminous than i/ersc/ieZ-detccted bump, star-forming 
DOGs, which have Lm = (3.1 ±0.4) x lO^^Lfr,, which i s 
consistent with the findings of IMelbourne et ahl (j2012[ ). 
Selection effects are investigated by calculating the IR lu- 
minosity of a representative CEOl template, scaled such 
that S24 = 100 ^Jy or S'250 = 8 mJy, as shown in Fig. [5l 
The local maxima in IR luminosity at z ^ 1.5 in the 
24 /im curve is due to the rest-frame 9.7 /im silicate ab- 
sorption feature. The lack of DOGs below the 24 ^m 
and 250 ^m limit at a given redshift leads us to conclude 
that the apparent trend in Fig. [5] that IR luminosity is 
increasing with redshift is a selection effect. 

Figure |6] shows dust temperatures as a function of IR 
luminosities for DOGs, color-coded by redshift. The av- 
erage characteristic dust temperature is T^ust = (34 ± 
7) K and (37 ±5) K for 77ersc/ie^detected and undetected 
DOGs, respectively, /fersc/ie^detected power-law DOGs 
and bump DOGs have average Tdust = (37 ± 6) K and 
(35 ±7) K, respectively, and is consistent with each other. 
We consider sample selection effects in the Tdust — ^ir 
plane by calculating IR luminosities and dust temper- 
atures of the CEOl templates with S'250 = 8 mJy and 
5*24 = 100 ^Jy. The 250 /im flux limit causes two se- 
lection effects: the first biases against hot sources, and 
the second is a lower luminosity limit that is redshift de- 
pendent (Fig. [6]). This second effect causes an apparent 
trend that the warmest and IR luminous sources lie at 
high redshifts. The 24 /im flux limit produces a similar 
effect. We observe that the dearth of sources at high lu- 
minosities and low dust temperatures is not caused by 
selection bias and i s a re al effect. This is consistent with 
iSvmeonidis et al.l (I2013D. who found that col d cirrus- 
dominated SEDs (jRowan-Robinson et al.l[2010[ ) are rare 
in the most IR luminous galaxies. 

The results shown in Fig. \6\ suggest that DOGs span 
a wider range of dust temperatures than z r^ 2 galaxy 



The errors are from bootstrapping. 
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Figure 5. IR luminosity as a function of redshift for Herschel- 
dctccted DOGs and the median IR luminosity for stack ed DOGs. 
A representative template from ICharv fc Elbad l|2001l ) scaled to 
the DOG 24 /im (solid curve) and 250 /xm detection limit (dashed 
curve) are also shown. A typical error bar is shown at the top right. 
The apparent trend that IR luminosity increases with redshift is a 
selection effect. 



Table 2 

Average DOG IR Luminosities and Dust Temperatures 



Type 


im 


idust 




(xIQI^Lq) 


(K) 


Detected^ 


2.8 ±0.3 


33 ±7 


Undetected" 


0.6 ±0.1 


37 ±5 


Bump* 


4.5 ±0.4 


34 ±7 


Power-law* 


3.1 ±0.4 


37 ±6 



Notes 

" Detected sources satisfy: S'250 > 3(t; undetected sources have 
•S'250 < 3(7. 

* Power-law (AGN-dominated) and bump (star-forming) median 
IR luminosities are from the detected sample. 

populations selected at longer wavelengths. While this 
may be true, it is more likely due to the different selec- 
tion effects associated with each galaxy population. For 
example, SMGs (by which we mean 850 /im or 1 mm se- 
lected sources) are biased towards detecting cold-dust 
dominated sourc es because hot sources are missed b y 
sub-mm surveys (jCasev et al.ll2009t iChapin et al.ll2009D. 
Optically-faint radio gala xies (OFRGs; iCasev et al.l 
l2009t iMagnelli et al.l |2010|) , which share similar stellar 
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Figure 6. Dust temperatures and infrared lumino sitie s for DOGs 
comp ared to o ther z ^ 2 ga l axy po p ulations (SMGsilMa gnclli ct a]] 
[20T^ . O FRGs:ICaseY et al.l 12^091 : IMagnelli et"an 12010.. 25 ^lra 
sources: [Casev et al.ll20Tir iRAC peakers: IMagdis et al.ll20Tol) . We 
show ffersc/iei-detected DOGs, colored by redshift. The average IR 
luminosity and dust temperature of i?ersc/ieZ-undetected DOGs at 
z = 2\s represented by the solid cyan circle. The black curves rep- 
resent the selection effects from estimating the d ust temperature 
and I R luminosity by considering templates from lCharv fc Elbad 
11200 II ) at 2 = 1.5, 2.0, and 2.5, with S250 = 8mjy, our 3cr detec- 
tion limit. The red curve is for S24 = 100 /ijy at 2 = 1.5. We 
conclude that the apparent trend that hotter sources are at higher 
rcdshifts may be a selection effect and that the absence of sources 
with high luminosities and colder temperatures is not a selection 
effect. 



masses, radio luminosities, and UV spectra as SMGs 
demonstrate this, while we also note that the radio- 
det ection limit is biased against the coldest sources 
(e.g lWardlow et ahl 120111) . DOGs are more insensitive 
to these selection biase s and thus show a wider range of 
temperatures. [K'lagdis et al.l (j2010t) found similar results 
when investigating the characteristic dust temperatures 
for IRAC peakers at z ~ 2 and showed that mid-IR se- 
lected sources bridge the gap in temperature ranges be- 
tween OFRGs and SMGs. We note that the 250 /.tm se- 
lected sources suffer from the same selection biases as our 
77ersc/ieZ-detected DOGs but shifted to higher luminosi- 
ties due to their shallower 250 ^va detection limit. 

3.2. Infrared Luminosity Function at z ^^ 2 

We compute the I R luminosity fu nction of DOGs using 
the 1/Kiax method (jSchmidtJUQTl) . defined as 



$(L)AL = 5] 



K. 



(4) 



where Knax is the maximum comoving volume of the ith 
source such that it would be detected and included in 
the sample. We consider the peak of the redshift dis- 
tribution using only DOGs at z = 1.5 — 2.5. For the 
Herschel-detected DOGs we use two flux limits to deter- 
mine T^nax: 5'24 = 100 /j,Jy; and 6*250 = 8mJy. These are 
the two detection limits of the survey. For the Herschel- 
undetected sample, the 24 /im flux limit alone was used 
to calculate Knax- The uncertainties are from Poisson 
statistics and binning errors, where the binning errors 
are calculated by generating the IR luminosity function 
1000 times from IR luminosities calculated from artificial 
SPIRE flux densities described in Section 13.11 and tak- 
ing the standard deviation per IR luminosity bin. The 



DOG IR luminosity function at z ^ 2 is presented in 
Fig. [7] and Table |3l The faint end of the IR luminos- 
ity function for Herschel-undetected and detected DOGs 
are coadded, which affects the lowest luminosity bin for 
Herschel-detected DOGs the most, showing a 0.20 dex 
increase. 

For comparison, the DOG IR luminosity function for 
77ersc/ie/-detected DOGs and all DOGs, calculated by ex- 
trapolating the infrared luminosity from 5*24 using CEOl 
templates, is also shown in Fig. [7l The contribution to 
the number density of /fersc/ie/-detected DOGs relative 
to the whole DOG population is never 100% at each IR 
luminosity bin, which means that some IR luminosities 
of the undetected population are severely being overes- 
timated. We find that the IR luminosities using this 
method are overestimated by a median factor of 1.8. This 
is consistent with t he previous studies of 24 /xm-selected 
galaxies at z ~ 2 (iHouck et al.l | 2005[: lYan et all 1 2007 ; 



baddi et al.l[200^ iPapovich et al.ll2007tlPope et al.ll2008 : 
Elbaz et al.ll2011l) and stresses the importance of far-IR 
observations, further highlighting that the use of 24 jim 
flux density to estimate IR luminosity is problematic at 
z > 1, due to features from PAH emission, silicate ab- 
sorption, and a power-law component from AGN. 

We compare the number densities of DOGs to the par- 
ent population of sources with 6*24 > 100 /iJy (Fig. [7]). 
The luminosity function of 24 /im sources is calcu- 
lated self-consistently using SPIRE data. There are 5932 
sources in COSMOS with ^24 > 100/iJy at z = 1.5-2.5, 
of which 32% are DOGs. Fig. [7] shows that DOGs have 
a smaller overall normalization in their luminosity func- 
tion (since they are fewer in number) and their relative 
contribution to the 24 /zm number density is roughly con- 
stant with increasing luminosity. Our results agree with 
the findings of Pcnner et al. (2012), where it was found 
that the distribution of IR luminosities for DOGs and 
the parent population of 24 /im sources are statistically 
indistinguishable . 

Figure [7] also sh ows the z = 1 — 3 SMC luminosity 
function from Wa rdlow et al.l ()2011l ) from two redshift 
bins. DOGs are more common than SMGs at z ~ 2, 
which is reflected in the higher normalization in the lu- 
minosity function, although in the HLIRG regime, SMGs 
dominate. This is consistent with the picture that in 
order for DOGs to represent an evolutionary stage af- 
ter the peak of star-formation rate in the SMC phase 
(jNaravanan et a l. 2010), they must have lower IR lumi- 
nosities and star-formation rates. As is shown in fig. |8l 
the total star-formation rate density (psfr) provided by 
the two populations are approximately even despite the 
number and intensity of sources. 

To calculate the contribution of DOGs to the psFR of 
the Universe at z ^ 2, we integrate the IR luminosity 
function and use equation |3l Figure [8] shows DOGs com- 
pared to other z ~ 2 galaxy populations. The total un- 
certainty in psFR is calculated from the quadrature sum 
of individual star-formation rate uncertainties and the 
standard deviation of psFR from the mock catalogs dis- 
cussed in Section r3. II Horizontal error bars represent the 
considered redshift interval. The value of psFR for DOGs 
at z = 1.5 - 2.5 is (3.2 ± 0.5) x lO'^ Mgyr-^Mpc"^ 
which contributes to 12 — 2 9% of the overall p s fr sA . 
z = 2 based on the models of iHopkins fc BeacomI (|2006f ) 
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Figure 7. IR luminosity function of DOGs in the COSMOS field at z ~ 2. Individually Herschel-detected DOGs and the results from 
stacking undetected DOGs are shown. We com pare this to an IR lumi nosity function fo r Herschel-detected D OGs and all DOGs generated 
from 24 ^m extrapolation using templates from lCharv fc Elbad {200j), classical SMGs l|Wardlow et al.ll20111 ). and 24 /im selected galaxies 
with 524 > 100 fi3y. The results of stacking allow us to estimate the faint end of the LF and we note that using 24 /ivci fiux density to 
calculate IR luminosity results in overestimation. DOGs have a higher normalization, "IJ*, but a lower luminosity turnover, L*, than SMGs. 



Table 3 

The IR Luminosity Function for DOGs at z ■ 



1.5-2.5 



logloC-t'lR/i©) 



<S> 



N" 



(Mpc~^dex~^) 



11.40- 11.80 
11.80- 12.15 
12.15- 12.50 
12.50 - 12.85 
12.85- 13.20 



(-3.79 ± 0.02)^ (660)'' 

-4.17 ± 0.21 (-3.98 ± 0.22)'' 150 (433)'' 

-4.03 ± 0.03 (-4.02 ± 0.09)'' 522 (522)'' 

-4.49 ± 0.03 252 

-5.40 ±0.09 31 



Notes 

"• N is the number of sources per luminosity bin. 
*" Numerical values in parenthesis include the stacked contribu- 
tion. 



and IBurgarella eTaLl (poTl . While 32% of the 24 ^m 
population at ^ = 1.5 — 2.5 arc DOGs, their fractional 
contribution to the star-formation rate density is slightly 
larger because of their typically higher IR luminosities. 
When comparing against all z = 1.5 — 2.5 sources with 
524 > 100 ^Jy, DOGs contribute 37% to the 24 //m psFR- 
This is consistent with the relative distribution of DOGs' 
infrared luminosities compared against the 24 /im popu- 
lation observed in the IR luminosity function. 

The undetected and power-law sources provide non- 
dominant contributions to the total psFR of DOGs. The 
77ersc/ieZ-undctected DOGs contribute 18% and power- 
law DOGs contribute just 9%. We note that even though 



power-law DOGs arc thought to be dominated by AGN 
emission in the IRAC bands, their far -IR emission is still 
likely dominated by star- formation (|Elbaz et al.l[2010i) . 
Indeed, even studies of the most active AGN have re- 
vealed that SED fits for ifersc/ie/-detected AGNs always 
required a star burst component in order to appear bright 
in the far-IR (jHatziminaoglou et al.l l2010( ) . As an at- 
tempt to quantify this claim, we use a simplified method 
to calculate an upper limit on the AGN contribution to 
the IR luminosity and star-formation rate in power-law 
DOGs and hence the contamination of asfb by AGN. We 
begin by scaling the AGN SEDs from iKirkpatrick et alJ 
(|2012[ ) to the 24 /im flux density of each power-law DOG 
and calculate the luminosity from the warm dust compo- 
nent. Then, by assuming that the warm dust component 
is entirely AGN-dominatcd and the cold dust component 
is entirely star-formation dominated, we can subtract the 
warm IR luminosity from the CEOl IR luminosity to cal- 
culate the residual contribution from star-formation. We 
find that power-law DOGs each have a maximum av- 
erage contribution of 70% to the IR luminosity, which 
could contaminate psFR by '^ 0.2 x 10~^ M0yr~^Mpc~^, 
which is only 6% of the total DOG psFR- In addition, 
we also estimate the dispersion of AGN c ontribution 
by no r malizing quasar SED t emplates fromlElvis et al.l 
(pOl . lRichards et al.l (J20M ). |PoIIe"tta et al.l (|2007| ) and 
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Figure 8. Star formation rate density (psfr) of all DOGs, Herschel-detected DOGs (offset by ^ = —0.1, for clarity), bump DOGs (offset 
by 2 = +0.1, for clarity) and all 524 > 100 /jjy sources in COSMOS at z = 1.5 — 2.5. We also show psfr for bump DOGs (524 > 100 
pjy) at 2 = 2 in the GOODS field from [Pope et al.l 12001) arid SMG s from lWardlow et al.l ll20Tll ) The evolution of psFR as a function of 
redshift from lHopkins &: Beacoriil 1 )20061 ) and IBurgarella et al.l 1 )20131) are also shown. Based on these models, DOGs contribute 12 — 29% 
to the total psFR of the Universe at 2 ~ 2. 



IDai et al] (|2012[ ) to the average power-law DOG 24 /^m 
flux density at z = 1.5 — 2.5 and assume that the SEDs 
have no emission associated with star-formation. Un- 
der this assumption, the average AGN contributions to 
the individual galaxies' IR luminosities range from 5% 
to 65%, depending on the SED, which corresponds to 
0.005% to 6% contribution to the total DOG psFR- 

We note that Pope et al. (200^ also examined bump 
(star-forming) DOGs at z '-^ 2 and they calculated 
Psfr ~ 1 x 10~^ Moyr~-^Mpc^'^. This value is lower 
than the bump psFR = 1-9 ± 0.3 x 10"^ Meyr^^Mpc"^ 
that we measure. However, it is difficult to determine 
wheth er these two values are significantly different be- 
cause IPope et all (|2008l ) do not provide an error on 
their measurement. We use their reported fractional 
error on the average IR luminosity (1.1 ± 0.7 x 10^^) 
to estimate that the minimum error on their psfr is 
~ 0.6 X 10^^ M0yr^^Mpc^^, in addition to the contribu- 
tion from the counting error from their 62 sources (com- 
pared to our 1137 bump sources at z — 1.5 — 2.5). We 
also note that the selection criteria for the two studies 
are sli ghtly different and if we were to use the IPope et "ahl 
(|2008|), we would identify 100 fewer bump DOGs ( 9% 
of bump DOGs are at z = 1.5 — 2.5). Therefore, we 
conclude that the two results are consistent, with our 
measurement providing a more accurate determination. 



3.3. Stellar Mass Build-up 

Using the stellar masses derived in lllbert et al.l (|20100 
(corrected to assume a Salpeter IMF by adding -1-0.24 
dex), and our derived star- formation rates using Her- 
schel data, we investigate where DOGs lie in the star- 
formation rate - stellar mass (SFR — M») plane. Disk 
galaxies with a steady star-formation mode are observed 
to form a tight correlation in their star-formation rates 
as a function of stellar mass, defining a "main sequence" 
(jDaddi et al.ll2007l:lElbaz et al.ll201lD . Outliers in this re- 
lati on are thought to be m erger-driven starburst galax- 
ies (jRodighiero et all 120111 and references therein). In 
the top panel of Fig. [9] we show the star- formation rates 
and stellar masses for //ersc/ieZ-detected DOGs, consid- 
ering only those at z = 1.5 — 2.5 to minimize the effects 
of redshift evolution. Average error bars are plotted for 
star-formation rates and the uncertainties in stellar mass 
are fixed to 0.5 dex, which covers the systematic offset 
range due to the choice of extinction laws and stellar 
population synthesis models. 

Fig. ini shows that power-law DOGs and bump DOGs 
cover the same ranges in stellar mass and star-formation 
rate in the SFR — M^ plane, as expected if the far-IR 
is star-formation dominated. Our findings are also con- 
sistent with previous studies that investigated the simi- 
larities in properties of far-IR SEDs of i/erscfteZ-selected 
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Figure 9. Star-formation rate (SFR; top panel) and specific star- 
formation rate (sSFR; bottom panel) as a function of stellar mass 
for DOGs a,t z = 1.5 — 2.5. Power-law DOGs and bump DOGs are 
statistically indistinguishable in the SFR — Af» plane. Th e sSFR at 
2 = 2 using the relation for star-forming galaxies from lElbaz et ahl 
(|2011h and its conversion to SFR for the displayed range of masses 
is shown as the thick solid orange line in both panels. The orange 
dotted lines represent a factor of two dispersion from the derived 
SFR and sSFR. DOGs have a large scatter about the main se- 
quence relation, having sources in the starburst, main sequence, 
and passive galaxy regimes. The thin horizontal black line in the 
top panel represents a minimum detectable star-formation rate at 
2 ~ 2, caused by the 24 fim flux density limit. Converting this to an 
sSFR value results in the diagonal line in the bottom panel, lead- 
ing us to conclude that the apparent negative correlation between 
sSFR and stellar mass is a selection effect. 



star-forming galaxies and AGN (iMullanev et al.ll2012f). 

The infrared main sequence from lElbaz et al.l (|2011[ ) 
for _ffersc/ieZ-selected star-forming galaxies at z = 2 is 
also shown in Fig. |9l DOGs have a significant amount 
of scatter about this relation, with 46% within a factor 
of 2 of the main-sequence, 24% above it and consistent 
with starbursts, and 31% below it in the more quiescent 
regime. 

The bottom panel of Fig. ^ shows the specific star- 
formation rate (sSFR = ^g^) as a function of stellar 
mass. The sSFR quantity measures the efficiency of star- 
formation, with its inverse giving a characteristic time- 
scale for the current burst length of star-formation ac- 
tivity. An apparent negative correlation in which lower 
mass DOGs exhibit higher sSFRs than their higher mass 
counterparts is observed for DOGs at 2; = 2, however 
this is largely a selection effect due to the fiux limit of 
our sample. On the top panel of Fig. [9] we use the mini- 



mum IR luminosity at z = 2 from our sample to represent 
a minimum detectable star-formation rate limit, shown 
as the horizontal line. We convert this to an sSFR value 
for a range of masses and this is shown as the diagonal 
line in the bottom panel. The logarithmic inverse age of 
the Universe in Gyr atz = 2is«— 9.5 (dashed line in 
Fig.[n]) and most DOGs have sSFRs larger than this, in- 
dicating that the observed phase of star-formation could 
be responsible for their observed stellar mass. 

Finally, we use the known redshift distribution and 
the sSFRs of DOGs to compare their volume densities 
to their likely progenitors, SMGs. We estimate that the 
volume density of observed DOGs at z = 1.5 — 2.5 is 
8 X 10^5 Mpc"^ Using the median DOG sSFR to as- 
sume a characteristic lifetime of the DOG phase to be 
^ 1 Gyr, we can correct this density for the burst duty 
cycle to derive a volume density for the progenitors to 
be ~ 1 X 10"'' Mpc"^. This is consistent to the S MG 
volume density derived from IWardlow et al.l (j2011[) at 
z = 1.5 — 3, which assumes the lifetime of the SMG 
phase to be 100 Myr, 10 times shorter than for DOGs. 
In this scenario, DOGs would have the sa me descendants 
as z ^ 2 S MGs, which, as dis cussed in IWardlow et al.l 
(|20T1 and iHickox eTall (|2012[ ). are likely to be 2-3L* 
early- type galaxies. 

4. CONCLUSIONS 

We use Herschd HerMES data in COSMOS to study 
the far-IR emission from DOGs. The main findings are: 

1. Out of 3077 DOGs, 51% are detected in Herschel 
{S250 > 3cr = 8 mJy) . We use stacking to probe the 
remaining undetected population and the stacked 
S'250 is on average a factor of ^ 1.5 fainter than the 
individual detection limit. 

2. The IR luminosity function of DOGs at z = 1.5 — 
2.5 is calculated and shows that they contribute 
30% to the luminosity function of 24 /im-selected 
galaxies. The stacked infrared luminosities provide 
significant contribution in the lowest detected IR 
luminosity bin, causing an increase of ^ 0.2 dex. 
IR luminosities derived from extrapolating 24 ^m 
flux densities of local galaxy templates are over- 
estimated by a factor of around 2 and therefore, 
direct observations in the far-IR that sample the 
cold-dust emission from galaxies should always be 
used, if available, to calculate more accurate esti- 
mates. 

3. DOGs contribute 10 — 30% to the overall star- 
formation rate density of the Universe and 39% for 
all 24 /.(m galaxies with 6*24 > 100 //Jy. We also 
note that when compared to the total DOG psFR, 
power-law (AGN dominated) DOGs provide minor 
contributions. The psFR for DOGs and SMGs are 
comparable at z ~ 2, however we note that DOGs 
are more numerous, with individually lower star- 
formation rates for DOGs than SMGs. 

4. DOGs have a large scatter in the SFR — M* plane, 
having sources in the starburst, main sequence 
and more quiescent galaxy regimes. The observed 
phase of star-formation for most DOGs is likely re- 
sponsible for their observed stellar mass. 
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